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We numerically investigate elementary excitations of the Heisenberg alternating-spin chains with 
two kinds of spins 1 and 1/2 antiferromagnetically coupled to each other. Employing a recently 
developed efficient Monte Carlo technique as well as an exact diagonalization method, we verify the 
spin-wave argument that the model exhibits two distinct excitations from the ground state which 
are gapless and gapped. The gapless branch shows a quadratic dispersion in the small-momentum 
region, which is of ferromagnetic type. With the intention of elucidating the physical mechanism 
of both excitations, we make a perturbation approach from the decoupled-dimer limit. The gapless 
branch is directly related to spin l's, while the gapped branch originates from cooperation of the 
two kinds of spins. 

PACS numbers: 75.10.Jm, 05.30.-d, 75.40.Mg, 75.30.Ds 



I. INTRODUCTION 

Extensive efforts have so far been devoted to verifying 
Haldane's conjecture H that the one-dimensional spin- 
s' Heisenberg antifcrromagnet exhibits qualitatively dif- 
ferent properties according to whether S is integer or 
half odd integer. The nontrivial energy gap immediately 
above the ground state was precisely estimated using 
numbers of numerical tools not only in the spin-1 case 
but also in the spin-2 case |5||| , while the valence- 
bond-solid model |tJ introduced by Affleck, Kennedy, 
Lieb, and Tasaki significantly contributed to the un- 
derstanding of the physical mechanism of the so-called 
Haldane massive phase. On the other hand, develop- 
ing the 0(3) nonlincar-er-model quantum field theory 
0, Affleck H pointed out that even integer-spin chains 
should be critical if a certain interaction is added to the 
pure Heisenberg Hamiltonian. Actually various numeri- 
cal methods p-|l5| revealed that the spin quantum num- 
ber is no more the criterion for the critical behavior in 
a wider Hamiltonian space. Recently several authors 
p6| , p7t even suggested the appearance of the Haldane- 
gap phases in half-odd-intcgcr-spin chains with a mag- 
netic field applied. Thus the low-temperature properties 
of one-dimensional quantum antiferromagnets with one 
kind of spins have more and more been elucidated. 

In such circumstances, there has appeared brand-new 
attempts jl8| |30| to explore the quantum behavior of 
mixed-spin chains with two kinds of spins. These studies 
are further classified according to their main interests. 
Several authors |l8|-pl| have been devoting their efforts 
to finding quantum integrable Hamiltonians and clarify- 
ing their critical behavior. Although the models consid- 
ered are generally complicated, the generic description of 
a certain family of Hamiltonians is interesting in itself 



and even allows us to guess the essential consequences 
of mixed-spin chains. A distinct attention is directed 
to mixed-spin chains with the simplest interaction be- 
tween the two kinds of spins. Alternating-spin Heisen- 
berg antiferromagnets with a singlet ground state [^6 29 
again present us the nontrivial gap problem |^,^]. Re- 
cently Fukui and Kawakami p6| have made a nonlinear- 
cr-model approach for a few models of this kind and have 
discussed a generic criterion for the critical mixed-spin 
chains. Their finding may stimulate many theoreticians 
to numerically investigate a variety of mixed-spin Hamil- 
tonians and even lead to synthesis of novel mixed-spin- 
chain materials. On the other hand, considering that all 
the mixed-spin-chain compounds synthesized so far ex- 
hibit a finite ground-state magnetization |3l[], we take 
a great interest in alternating-spin Heisenberg antiferro- 
magnets with ferrimagnetic ground states. This is the 
subject we discuss in the present article. 

Let us introduce a Hamiltonian of alternatively aligned 
two kinds of spins S and s which are antiferromagneti- 
cally coupled to each other: 



n 



N 



Sj + fi.s, ■ Sj+i) , 



(1.1) 



where S 2 = 5(5- 



1), s 2 = s(s+l), and N is the number 



of unit cells. The bond alternation 5 has been introduced 
for a discussion presented afterwards. We assume 5 > s 
in the following without losing generality. Because of the 
non-compensating sublattice magnetizations, this system 
exhibits the ferrimagnetism instead of the antiferromag- 
netism. Appl ying the Lieb-Mattis theorem |32| to the 
Hamiltonian (1.1), we immediately find (5 — s)iV-fold 
degenerate ground states. The Goldstone theorem |3^] 
further allows us to expect a gapless excitation from the 
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ferrimagnetic ground state. Therefore we here take lit- 
tle interest in the simple problem whether the system is 
gapped or gapless. Alcaraz and Malvezzi |||] investigated 
the two cases of (S, s) = (1, 1/2) and (S, s) = (3/2, 1/2) 
and actually showed that in both cases the chain is de- 
scribed in terms of the c = 1 Gaussian conformal field 
theory. Suggesting that this should be the generic sce- 
nario for arbitrary ferrimagnetic Heisenberg chains, they 
further predicted the appearance of quadratic dispersion 
relations at the Heisenberg points, which means that the 
model possesses a ferromagnetic character. On the other 
hand, applying the spin-wave theory to the model, two 
groups [ 24p5|] have recently predicted that there exists a 
gapped branch of elementary excitations as well as a gap- 
less ferromagnetic branch. Their prediction is quite inter- 
esting because it implies the coexistence of the ferromag- 
netism and the antiferromagnetism in the ferrimagnets. 
This is the motivation for the present study. Employing 
a quantum Monte Carlo technique and an exact diago- 
nalization method, we here calculate energy eigenvalues 
of the elementary excitations with the total magnetiza- 
tion Y,f=i(Sj + s f) = M = ( S - S ) N =F 1, which cor- 
respond to the ferromagnetic and the antiferromagnetic 
branches, respectively. The numerical results are com- 
pared not only with the spin-wave calculation but also 
with a perturbation approach from the decoupled-dimcr 
limit (5 = 0). 



II. THE SPIN- WAVE APPROACH 



First, we briefly review the spin- wave-theory result 
p4| , p5| , which allows us to have a qualitative view of the 
low-energy structure. We start from a Neel state with 
M = (S — s)N, namely, we define the bosonic operators 
for the spin deviation in each sublattice as 



E = -J(l + S)SsNJ 



k 



Sj~ = V2Sa,j, Sj = \/2S oj , Sj = S — a]aj , , , 
F,t c--—./oZh. ,2- _ D _i_f,ti. I ■ > 



2sb), sj 



'2s bj , s z j 



S + bjbj 



In order to obtain the dispersion relations of the spin- 
wave excitations, we handle the boson Hamiltonian up 
to quadratic order. We define the momentum represen- 
tation of the bosonic operators as 



2iajfe t 



N 
N 

; t _ 1 2ia7fc/,t 
' - ^Z^ C °3 

3 = 1 



(2.2) 



where k = irl/Na (/ = -N/2+1, -N/2+2, N/2) with 
a being the distance between two neighboring spins. We 
note that here the unit cell is of length 2a. Carrying out 
a Bogoliubov transformation 



Oik 



-iaX k /2 cosh g k Qk 



e iaA fc /2 sinh6)fc 



(3 k = e iaAfc / 2 sinh# fe 4 + e- iaAfc / 2 cosh6» fc b k . 



(2.3) 



with 



(1 + S) tan[(& - X k )a] -(1-6) tan(afc) = , (2.4) 



2VSs 



tan(2 * fe) a+s)(s+ S ) 

x \J(1 + S) 2 cos 2 (ak) + (1 - S) 2 sin 2 (ak) , (2.5) 
we reach the diagonal Hamiltonian, 

H = E + J2 ("k a l a k + . ( 2 - 6 ) 

fe 

where 



165 Ss sin 2 (ak) - (1 + 6)(S + s) 



(2.7) 



y/(l + S) 2 (S - s) 2 + 165Ss sin 2 (ak) + (1 + S)(S — s) 



(2.8) 



Thus the spin- wave approach suggests the coexistence of 
the ferromagnetism and the antiferromagnetism in the 
present system. In the ferromagnetic branch (u> k ), the 
spin wave reduces the total magnetization and exhibits 
a quadratic dispersion relation in the small-momentum 
region: 



fe — o 



JSSs(2ak) 2 

(1 + *)(£-«)' 



(2.9) 



while in the antiferromagnetic branch (oJ k ), the spin 
wave enhances the total magnetization and have the 
gapped excitation spectrum: 



J(l + 6)(S-s) + 



J5Ss(2ak) 2 
(l + 5)(S-s) 



(2.10) 



We note that the qualitative character of the model re- 
mains unchanged over the whole region of 5. It is the 
purpose in this article to confirm this scenario employing 
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efficient numerical methods. 



III. A PERTURBATION APPROACH 

Second, prior to the numerical approach, let us carry 
out a perturbation calculation with the intention of elu- 
cidating the nature of the elementary excitations. In the 
decoupled-dimer limit, we can easily find the low-lying 
eigenstates. Figures 0(a), 0(b), and |l](c) represent, re- 
spectively, the ground state, \^f), the ferromagnetic exci- 
tation at an arbitrary unit cell j, l^j), and the antifer- 
romagnetic excitation at an arbitrary unit cell j, I'J'j - ), 



of the Hamiltonian (|L1|) with (S, s) = (1, 1/2) at 5 = 0. 
When we turn on the exchange interaction between the 
dimers, the localized excitations can hop to neighboring 
unit cells with an amplitude proportional to S. We take 
account of this effect using the degenerate perturbation. 
We introduce a representation of matrix-product type as 



,1 



' 9 N S ® h% . 



I*/) = gl s ® K 



where 



(3.1) 

® 9n ® h N , (3-2) 
■®9N®h s Nl (3.3) 



both branches remain unchanged as 5 increases: The fer- 
romagnetic branch tu^ is gapless and proportional to k 2 
in the small-/c region, while the antiferromagnetic branch 
uj~1 is gapped. 

For an arbitrary combination of (S, s), the similar ar- 
gument can be developed and qualitatively the same re- 
sult is obtained. For example, in the case of (S, s) — 
(3/2, 1/2), we find the dispersion relations 

uj- = ljS[l-cos(2ak)] + 0(S 2 ), (3.8) 
8 

oj+ = 2J+ -JS[7 - 6 cos(2afc)] + 0(S 2 ) , (3.9) 
8 

with Eq — — (5/4) J(l + S/A)N, considering the elemen- 
tary excitations shown in Fig. |^. Thus we are more and 
more convinced that the scenario of the low-energy struc- 
ture should be valid for an arbitrary Heisenberg ferrimag- 
net. Actually, recent numerical studies [§2],|4} on the 
present Hamiltonian reported that the low-temperature 
properties are essentially the same regardless of the val- 
ues of S and s as long as they differ from each other. 
Alcaraz and Malvezzi found that the model with 
exchange anisotropy is described in terms of the Gaus- 
sian critical theory in both cases of (S, s) — (1, 1/2) and 
(5, s) = (3/2, 1/2). In such circumstances, we restrict our 
numerical investigation to the case of (S, s) = (1, 1/2). 



9j = 



[10), V2\+)j 
[-V2~\+)j 
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the s 



IV. NUMERICAL PROCEDURE 

In order to calculate the low-lying eigenvalues of the 
model, we here use two numerical tools, which possess, 
respectively, both advantageous and weak points of their 
own and are complementary to each other. 

One is the exact-diagonalization method employing the 
Lanczos algorithm. Calculation of the energy levels re- 
±)j, \0)j being the SJ-eigenstates and | t)j, I i)j duces to the diagonalization of the 6 N x 6^ matrix rep- 



"-IT)/ 




"-It)i" 


. 1 ih . 








(3.4) 



(3.5) 



-eigenstates. Now the dispersion relations of the 



eigenstates with M = N/2 =p 1, oJf , are calculated as 



E G 



-J5[l-cos(2ak)} + 0(5 2 ) 



(3.6) 



3 1 cr 
2 J+ 18 J ^ 7 



-Eg 

12cos(2a/c)] + 0(<5 2 ), (3.7) 



where 
N 



I* 



kl 

1/2 V^A' 
2^3 = 1 



= n-v 2 y: 



J = l 

are 



1^) 



forms of l^j) and 



and 

the Fourier trans- 

respectively, and Eg = 
*|W|*)/(*|*) = -J(l + 5/9)N is the ground-state 
energy within the first order of S. Although the Heisen- 
berg point (6 = 1) is far from the decoupled-dimer limit, 
it is interesting enough that the qualitative characters of 



resenting the Hamiltonian. In constructing basic states, 
we use direct products of the single-spin states indicated 
by the projection values. The present Hamiltonian com- 
mutes with the total magnetization M and therefore 
splits into 3iV + 1 blocks labeled by M (on the assump- 
tion that N is even). Since we perform the calculation 
under the periodic boundary condition, each eigenvalue 
is further classified by its total wave number k. 

We treat the chains ofiV = 8,10,12, where we restrict 
the calculation to the lowest energy level in each sub- 
space. Since our main interest is to reveal the nature of 
the elementary excitations, we investigate the subspaces 
of M = N/2 and M = N/2 =F 1 for all the chain lengths 
we treat. We calculate the subspaces of M < N/2 — 2 
as well for N = 8, 10 in order to further elucidate the 
ferromagnetic nature of the model. Due to the two kinds 
of spins, the construction of the basic states is somewhat 
complicated. The base dimension reaches ten million in 
the case of N = 12 and M = N/2 - 1. 

The chain length we can reach with the diagonaliza- 
tion method is much smaller than one with a Monte Carlo 



3 



technique. However, having in mind the small correlation 
length of the present system |24|,^5|], the diagonalization 
result is fruitful enough to obtain a general view of the 
low-energy structure. Actually, the ground-state energy 
for N — 12 coincides with the N — > oo extrapolated 
value within the first five digits. We further note that 
even the ground-state energy for N = 8 is so close to the 
thermodynamic-limit value as to show the coincidence 
within the first three digits. This fact is quite helpful 
in estimating the dispersion relations in the long-chain 
limit, although the eigenvalues for different chain lengths 
should be distinguished. 

The other approach is based on a quantum Monte 
Carlo technique |l2||34|-^6| which one of the authors 
has recently developed. The idea is summarized as ex- 
tracting the lower edge of the excitation spectrum from 
imaginary-time quantum Monte Carlo data at a low 
enough temperature. The imaginary-time correlation 
function S(q, r) is generally defined as 



S(q,r) = (c 



Ht 



0~e 



-Ht 



o z - q ) 



(4.1) 



where O q — N 1 2j=i Oje lq ^ a ° is the Fourier transform 
of an arbitrary local operator Oj with ao being the length 
of the unit cell, and (A) = Tr[e-P n A]/Tr[e-P n ] denotes 
the canonical average at a given temperature /3 _1 = k^T . 
While S(q, r) as a function of r generally exhibits a com- 
plicated multi-exponential decay, it may efficiently be 
evaluated at a sufficiently low temperature as jj4| 

S(q, t) = V 1(1; ko\S'\l; k + q) | 2 e -^(*°+«)-MM] , 



l 



(4.2) 



where \l;k) (I = 1,2,- ■■) and Ei(k) {E x (k) < E 2 (k) < 
• • •) are the Ith eigenvector and eigenvalue of the Hamil- 
tonian in the fc-momentum space, and fco is the mo- 
mentum at which the lowest-energy state in the sub- 
space is located. Now it is reasonable to approximate 
.Ei (fco + <?) — -Ei(fco) by the slope — dhi[S(q, t)]/5t in the 
large-r region satisfying 



-[E 2 (ko +q)- -Bi (fco + ?)] > In 



|(l;fc |5 g z |n;fc +g)| 2 
|(l;fc |5||l;fc + g)| 2 ' 

(4.3) 



for an arbitrary n. 

In case the lower edge of the spectrum is separated 
from the upper bands or continuum by a finite gap, or 
in case of its spectral weights \ k \O^\l; k + q}\ 2 being 



relatively large, the inequality (4.3) is well justified even 
at small r's and a logarithmic plot of S(q, r) is expected 
to exhibit fine linearity in a wide region of r. Actu- 
ally, for single-spin Heisenberg chains with an arbitrary 
spin quantum number or a certain bond alternation, it 
was shown ]T^j3^ -|3^] that the r-dependence of S(q,r) 
is essentially approximated by a single exponent at each 
momentum q at a sufficiently low temperature. 



Here, due to the two kinds of spins in a chain, Oj is not 
uniquely defined. We show in Fig. [| logarithmic plots 
of S(q, t) as a function of r setting several operators for 
Oj. We have set (/3J)" 1 and n equal to 0.02 and 200, 
which are, respectively, low and large enough |?6| to re- 
move the finite-temperature effect and the finite- n effect. 
In order to estimate each S(q,r), we have carried out a 
few million Monte Carlo steps spending several days on 
a supercomputer or a few weeks on a fast workstation. 
Energy difference between the ground state and the low- 
est state with an arbitrary q is obtained through S(q, r) 
calculated in the subspace of M = 0. 

In Fig. ^ we plot the excitation energies as a function 
of q obtained by estimating the slope — dhi[S(q, t)]/8t 
in the largest-r region available. In the case of Oj = s*, 
the multi-exponential behavior of S(q, r) is remarkable 
even in the large-r region and thus prevents us from pre- 
cisely estimating the energy eigenvalues. In all the other 
cases, we obtain useful data within the numerical preci- 
sion. With the present data, we can at least conclude that 
spin-l/2's do not have much effect on the lowest-lying ex- 
citations, which suggests that the elementary excitations 
of ferromagnetic nature may qualitatively be described 
by the simple picture shown in Fig. |l|(b) even at the 
Heisenberg point 5=1. While Oj = Sj+s* brings some- 
what higher energies than Oj = S* and Oj = — s|, 
we find no difference beyond the numerical accuracy be- 



tween the cases of O, 



and Oj = Sj 



It will be 



shown in the next section that the thus-obtained lowest- 
lying energy eigenvalues are in good agreement with the 
exact-diagonalization result. This method is applicable 
to rather long chains but is not successful in obtaining the 
higher-lying eigenvalues except for the special fortunate 
cases |34| ]. Thus the diagonalization technique is nec- 
essary and useful for the antiferromagnetic branch with 
M > N/2 even though it is inferior to the Monte Carlo 
method in treating long chains. 



V. NUMERICAL RESULTS AND DISCUSSION 



We plot in Fig. |5| the quantum Monte Carlo and 
the exact-diagonalization calculations of the excitation 
energies as a function of momentum q for the chain 
without bond alternation, where the results of the spin- 
wave theory and the first-order perturbation from the 
decoupled-dimer limit are also shown. The lower band is 
the lower edge of the excitation spectrum and consists of 
the lowest-lying eigenvalues with M = N/2 — 1. It ex- 
hibits a quadratic dispersion at small q's as was expected. 
The upper band consists of the lowest-lying eigenvalues 
with M = N/2 + 1 and is separated from the ground 
state by a finite gap. It is the scenario predicted by the 
analytic approaches that we here observe. The quan- 
tum Monte Carlo finding is in good agreement with the 
diagonalization result. The diagonalization calculation 
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indicates that the chain- length dependence of the disper- 
sion relation is quite weak even in the vicinity of the zone 
boundary and the zone center, which is consistent with 
the extremely small correlation length J24|,|2^] . 

In order to perform a quantitative comparison between 
the numerical findings and the analytic results, let us con- 
sider the curvature of the dispersion, v, which is defined 
by 

uj- = v{2ak) 2 . (5.1) 

The quantum Monte Carlo method, the spin wave calcu- 
lation, and the perturbation approach, respectively, give 
v QMG/j = o.37(l), v sw /J = 1/2, and v^/J = 2/9. 
The spin-wave theory overestimates the true value, while 
the perturbation approach underestimates that. We have 
made an attempt to obtain another estimate of v using 
the exact-diagonalization result. Although the Lanczos 
algorithm results in much more precise raw data than 
the Monte Carlo technique, yet the attempt was not so 
successful as one using the Monte Carlo data because of 
the lack of data points. However, we have confirmed that 
the diagonalization estimate of v possibly coincides with 
V QMC wjthjn the numerical accuracy. We note that the 
spin-wave estimate v accords with one for the Heisen- 
berg ferromagnet of spin 1/2 in the unit of the unit-cell 
length being equal to unity. Furthermore the calculation 
of the lowest levels in the subspace of M = TV/ 2— 2 results 
in energy eigenvalues bellow the two-magnon continuum. 
The obtained state should be a two-magnon bound state 
|57j. All these facts again emphasize the ferromagnetic 
aspect of the present model. However, in contrast with 
the ferromagnet, the true value «Q MC is reduced from 
v sw due to the quantum effect. Thus the elementary ex- 
citations in the sector of M < N/2 can be regarded as 
spin waves modified by quantum fluctuations. 

For the antifcrromagnetic branch, on the other hand, 
both spin-wave analysis and perturbation calculation are 
not so successful as ones for the ferromagnetic branch. 
Considering that quantum effects are, in general, much 
more remarkable in antiferromagnets than in ferromag- 
nets, the most naive picture of the elementary excitations 
illustrated in Fig. [I] and Fig. ^ may have to be signifi- 
cantly modified for the antiferromagnetic branch. 

Finally in this section, we show in Fig. ^| the calcula- 
tion for the chains with bond alternation. Although the 
quantum Monte Carlo calculation is generally in good 
agreement with the diagonalization result, the agreement 
seems to be somewhat poorer in the small-<5 region. This 
is convincing keeping in mind that the decrease of 5 may 
cause the freezing of the spin configuration in Monte 
Carlo sampling and therefore a huge number of Monte 
Carlo steps are needed to refine the data accuracy in the 
small-<5 region. It is needless to say that the perturbation 
calculation is more justified in the small-i5 region. As the 
model approaches the decoupled-dimer limit, both fer- 
romagnetic and antiferromagnetic bands become flatter 
and approach and (3/2) J, respectively. The spin-wave 



result is fairly good in the ferromagnetic branch but rela- 
tively poor in the antifcrromagnetic branch. This is con- 
vincing considering that the spin wave correctly describes 
the low-lying excitations of the fcrromagnets, while it is 
valid at most qualitatively for the antiferromagnets. 



VI. SUMMARY 

We have investigated the low-energy structure of the 
ferrimagnetic alternating-spin chains with spin 1 and spin 
1/2. Motivated by the spin- wave analysis of the model, 
we have mainly calculated the eigenvalues with an arbi- 
trary momentum of the one-magnon states, namely, the 
lowest-lying eigenvalues in the subspaces of M = N/2^fl. 
The chain-length dependence of the dispersion relations 
is extremely weak, which is consistent with the consid- 
erably small correlation length of the system, £ < 2a 
^4j,^5|. The qualitative character of the model remains 
unchanged under the existence of bond alternation. The 
ferromagnetic branch is gapless and shows a quadratic 
dispersion in the small-momentum region, while the anti- 
ferromagnetic branch is separated from the ground state 
by a finite gap A. A/ J in the thermodynamic limit is es- 
timated to be 1.75914(1) at the Heisenberg point 5 = 1. 
We made an attempt to understand the low-lying exci- 
tations through the first-order perturbation calculation 
from the decoupled-dimer limit. The ferromagnetic exci- 
tations are more or less dominated by spin l's, whereas 
the mechanism of the antiferromagnetic excitations was 
less revealed. Quantum Monte Carlo snapshots |38| may 
help us to inquire further into the antiferromagnetic fluc- 
tuations. 

We have further carried out the diagonalization calcu- 
lation in the subspace of M = N/2 — 2 and found the 
two-magnon bound state bellow the continuum. This 
fact emphasizes the ferromagnetic aspect of the sector 
of M < N/2 in the ferrimagnet. We expect that the 
Heisenberg ferrimagnet behaves like a ferromagnet at low 
enough temperatures, while its antiferromagnetic aspect 
may appear at k^T > A. 
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FIG. 1. Schematic representations of the ground state (a), 
the elementary excitation in the subspace of M — N/2 — 1 (b), 
and the elementary excitation in the subspace of M = N/2+1 
(c) of the ferrimagnetic chain with spin 1 and spin 1/2 in the 
decoupled-dimer limit. The arrow (the bullet symbol) denotes 
a spin 1/2 with its fixed (unfixed) projection value. The solid 
(broken) segment is a singlet (triplet) pair. The circle repre- 
sents an operation of constructing a spin 1 by symmetrizing 
the two spin 1/2's inside. 



FIG. 2. Schematic representations of the ground state (a), 
the elementary excitation in the subspace of M — N/2 — 1 (b), 
and the elementary excitation in the subspace of M — N/2 + 1 
(c) of the ferrimagnetic chain with spin 3/2 and spin 1/2 in 
the decoupled-dimer limit. The notation is the same as one 
in Fig. 1 except for the circle representing an operation of 
constructing a spin 3/2 by symmetrizing the three spin 1/2's 
inside. 



FIG. 3. Logarithmic plots of S(q, r) versus the imaginary 
time r at several choices of Oj taking a few values of mo- 
mentum q for the Heisenberg ferrimagnetic chain of N — 32 
with (S,s) = (1,1/2) and 5 = 1: Q (2aq/n = 16/64), + 
(2aq/n = 48/64) with Oj = S"\ V {2aq/n = 16/64), x 
(2aq/n = 48/64) with Oj = sj; O (2aq/n = 16/64), □ 
(2aq/n = 48/64) with Oj = S- + sj; A (2aq/n = 16/64), 
O {2aq/n = 48/64) with Oj = S] - s). The numerical un- 
certainty is all within the size of the symbols. 



G 



FIG. 4. Quantum Monte Carlo estimates of excitation en- 
ergies as a function of q for the chain of N = 32. Here Oi 
□ , O, and x have, respectively, been obtained from S(q,T)'s 
with 3 = S|, Oj = sj, 3 = S| + s z j, and 3 = S| - s| at 
the subspace of M = 0. The error bars are attached to the 
data obtained with Oj = and Oj = S| + sj. The numerical 
uncertainty of the rest of the data is all within the size of the 
symbols. GS denotes the ground state. 

FIG. 5. Quantum Monte Carlo and exact-diagonalization 
calculations of the lowest energies as a function of momentum 
in the subspaces of M = N/2 =p 1 for the Heisenberg ferrimag- 
netic chain of N = 32 with spin 1 and spin 1/2. The numerical 
uncertainty is all within the size of the symbols. The results 
of the spin-wave theory and the first-order perturbation from 
the decoupled-dimer limit are also shown by solid and broken 
lines, respectively. GS denotes the ground state. 



FIG. 6. The lowest energies as a function of momentum in 
the subspaces of M = N/2-1 (a) and M = N/2+1 (b) for the 
bond-alternating Heisenberg ferrimagnetic chains of N — 32 
with spin 1 and spin 1/2. Here O (I), D (-), O (+), and 
A (x) represent the exact-diagonalization (quantum Monte 
Carlo) estimates at S = 0.8, 8 = 0.6, 8 = 0.4, and 8 = 0.2, 
respectively. The numerical uncertainty is all within the size 
of the symbols. The results of the spin-wave theory and the 
first-order perturbation from the decoupled-dimer limit are 
also shown by solid and broken lines, respectively, where the 
values of 8 are 0.8, 0.6, 0.4, and 0.2 from top to bottom except 
for the perturbation result for the antiferromagnetic branch 
in the small-g region with the values of 8 being 0.2, 0.4, 0.6, 
and 0.8 from top to bottom. GS denotes the ground state. 
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